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ABSTRACT 

We report on a survey devoted for the search of exo-planets around young and nearby stars 
carried out with NACO at the VLT. The detection limit for 28 among the best available targets 
vs. the angular separation from the star is presented. The non-detection of any planetary mass 
companion in our survey is used to derive, for the first time, the frequency of the upper limit 
of the projected separation planet-stars. In particular, we find that in 50% of cases, no 5Mj 
planet (or more massive) has been detected at projected separations larger than 14 AU and no 
10Mj planet (or more massive) has been detected at projected separations larger than 8.5 AU. In 
100% of cases, these values increase to 36 AU and 65 AU respectively. The excellent sensitivity 
reached by our study leads to a much lower upper limit of the projected planet-star separation 
compared with previous studies. For example, for the (3 Pictoris group, (~12 Myr), we did not 
detect any 10Mj planet at distances larger than 15 AU. A previous study carried out with 4 m 
class telescopes put an upper limit for 10Mj planets at ^60 AU. For our closest target (V2306 
Oph - d = 4.3 pc) it is shown that it would be possible to detect a 10Mj planet at a minimum 
projected separation from the star of 1 AU and a 5Mj planet at a minimum projected separation 
of 3.7 AU. Our results are discussed with respect to mechanisms explaining planet formation and 
migration and to forthcoming observational strategies and future planet finder observations from 
the ground. 

Subject headings: binaries: close — planetary system — stars:low-mass, brown dwarfs 



1. Introduction 

At the present time more than 100 exo-planets 
are known which have been mostly discovered by 
the indirect method of the variation of the radial 
velocity induced by planets orbiting around the 
central star 2 . It can be shown that young planets 



1 Based on observations collected at the European 
Southern Observatory, Chile. Program 70. C — 0777D, 
70.C - 0777B and 71. C - 0029A 
2 Geneva Chronology Catalog 2004 - 
http: / / www.obspm.fr / encycl / catalog.html 



can be detected with adaptive optics (AO) imag- 
ing techniques at 8-10 meter telescopes if they are 
sufficiently distant from their parent star. Since 
it is reasonable to assume a coevality between the 
parent star and the planet, we can retrieve from at- 
mospherical models (Burrows et al. 1997, Baraffe 
et al. 2003) that planets having a mass in the 
range 3-10 Mj and orbiting around young late- 
type stars (10-200 Myr) have a typical brightness 
contrast with respect to the parent star of the or- 
der of 10 2 -10 6 , corresponding to a magnitude dif- 
ference of AM of 5-15 mag. This means that a 
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planet can be detected at a few tens of AUs from 
parent stars having distances < 50 pc. First re- 
sults on the search for young and massive planets 
are reported by Neuhauser et al. (1997) using 4 m 
class telescopes and Macintosh et al. (2001) and 
Kcisler et al. (2003) using 10 m class telescopes. 
Recently (Chauvin et al. 2004), a companion of ~ 
5±2Mj was detected at around 50 AU from the 
central star 2M1207, a very late type star (M8) 
in TWA Hydrae (~ 10 Myr). This object is, at 
the present time, the most interesting giant planet 
candidate ever detected and further confirmation 
of the orbital motion is needed. 

From the point of view of the scientific moti- 
vations related to this kind of search it is cer- 
tainly an important challenge to image directly 
a planet around another star, since it can pro- 
vide a direct evidence for the existence of such 
a planet. At the present time only a few exo- 
planets have been detected with the photometric 
technique called 'transit' (HD209458b - Charbon- 
neau et al. 2000, OGLE-TR-113 (Konacki et al., 
2004, Bouchy et al. 2004), OGLE-TR-56 (Torres 
et al., 2004), OGLE-TR-132 (Bouchy et al. 2004), 
OGLE-TR-111 (Pont et al., 2004) and TrRes-1 
(Alonso et al., 2004)). Moreover, direct imaging 
is complementary to other detection methods like 
the radial velocity variation technique. The latter 
is sensitive to objects orbiting at close distances 
from the parent star (the radial velocity induced 
by the presence of a companion is indeed propor- 
tional to the inverse of the square root of the or- 
bital distance) while the AO imaging presented 
here would be sensitive to planets orbiting at much 
larger distances. 

According to statistical results retrieved by ra- 
dial velocity studies (Marcy et al. 2003), most of 
the known exo-plancts orbit at distances smaller 
than 1 AU and the semi-major axes of exo-planets 
found so far are not larger than ^6 AU. Since the 
radial velocity technique is heavily biased towards 
short separations, the derived major-axis distribu- 
tion is certainly not representative. Not much is 
known at the present about potential planets that 
could exist at distances larger than about 6 AU. 
We note that the formation of exo-planets is quite 
unlikely at distances of a few tens of AU. At these 
large distances the density of the planetesimals is 
low and the gravity from the central star is weak. 
This makes the process of planet formation ineffi- 



cient and slow (Thommes et al., 1999, 2002). How- 
ever, different models predict the transfer (jump) 
or migration of planets to distances larger than 5- 
10 AU from the central stars. Among these mod- 
els, we recall the gravitational scattering model 
which claims that in multi-planetary systems, due 
to gravitational interactions between the planets, 
it can happen that the lowest mass objects are 
ejected on hyperbolic trajectories and some plan- 
ets can move on stable orbits of a few tens of AU 
from the parent star. Such models were proposed 
by Weidenschilling & Marzari (1996) and Rasio 
& Ford (1996) and revised more recently by Pa- 
paloizou & Terqucm (2001) and Ford et al. (2001). 
An alternative mechanism that could induce an 
outward migration is the planet-disc interaction 
(Veras & Armitage, 2004). It was proven that 
stable orbits of a few tens of AU can be attained 
by a few Jupiter mass planet under the effect of 
torques induced by the disc (i.e at early stages of 
the star formation) if the system has photoevapo- 
ration wind. Numerical simulations (Veras & Ar- 
mitage, 2004) showed that, under the assumption 
of a typical solar-mass photoevaporation wind of 
the order of 1-5 x 10 -9 MQyr -1 , an outward migra- 
tion of a planet with few Jupiter masses can take 
place with a frequency of the order of 25 % (with 
respect to a sample of a few hundreds of simula- 
tions). It is, therefore, interesting to find out if 
young and massive planets are present up to a few 
tens of AU from the central star. 

For these reasons, we planned and performed an 
extended survey aiming to search for massive exo- 
planets orbiting around nearby and young stars 
using one of the best available instruments at the 
present time for high-contrast AO imaging ob- 
servations: NACO at the Very Large Telescope 
(VLT) (see Section 3). 

The plan of the paper is the following. In Sec- 
tion 2 the criteria used for the selection of the tar- 
gets are presented. In Section 3 we describe the 
observational strategy and we provide technical in- 
formation on the observations. In Section 4.1 we 
present our data processing, in Section 4.2 we dis- 
cuss an analysis of planet-like features identified in 
deep images and in Section 4.3 we give a detection 
limit estimate for the observed targets. In Section 
5 we discuss our results and finally, in Section 6 
the conclusions of our study are presented. 
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2. Target selection 

The properties of planets that we are searching 
for are: (a) the contrast AM between the planet 
and the central star should be 5-15 mag (b) and 
the planets should orbit at distances a > 1 AU 
from the central star. 

The planet intrinsic luminosity, retrieved by 
evolutionary models (Burrows et al. 1997, Baraffe 
et al. 2003), depends on the age of the system and 
the mass of the planet. The magnitude of a planet 
depends also on the wavelength used for the obser- 
vations. Due to the fact that planets strongly fade 
with age we selected stars as young as possible. 
Following the evolutionary models, we define our 
'selection' ranges as follows: age ~ 10-200 Myr, 
mass of the planets ~ 3-10 Mj and wavelengths 
in the H and K s bands. Obviously the nearer the 
target is from the observer, the smaller is the dis- 
tance from the central star that we can reach and 
probe. Therefore, we gave priority to nearby ob- 
jects. We note that most (~ 70 %) of our targets 
have distances < 32 pc and only 3 objects are lo- 
cated between 50 and 77 pc. In order to optimize 
the contrast conditions, we gave priority to late- 
type stars (spectral type K-M). Only one target 
(HD 155555) has a spectral type G5IV. Table 1 
and 2 summarizes the properties of the selected 
targets and the theoretical magnitude of planets 
that might orbit around them. The magnitude of 
the stars in the H and K s bands are retrieved from 
the magnitude in V band and the spectral type 
(Simbad Astronomical Database 3 ). Table 1 gives 
also the references from which the targets were se- 
lected. A relatively large fraction of targets belong 
to the TW Hydrae local associations, (3 Pictoris 
moving groups and Tucana/Horologium associa- 
tion (Zuckcrman et al. 2001a, 2001b; Song et al. 
2003, Barrado y Navascues et al. 1999, Barrado y 
Navascues 2001). The age of TW Hydrae is esti- 
mated to be ~ 10 Myr (Stauffer et al. 1995), the 
age of j3 Pictoris is ~ 12 Myr (Barrado y Navascues 
et al. 1999) and the Tucana/Horologium age is ~ 
30-40 Myr (Zuckcrman & Webb 2000). Several 
other targets, members of young stellar kinematic 
groups, were selected from Montes et al. (2001). 

We want to emphasize that the ages given in 
Table 1 are often approximate values, particularly 

3 http://simbad. u-strasbg.fr/ 



for those stars which are not members of associ- 
ations or moving groups. For sufficiently young 
low-mass stars (< 30-50 Myr) their location in 
theoretical Hcrtzsprung Russel diagrams can be 
used for an age estimate. For older stars (> 50 
Myr) the isochrones arc often indistinguishable 
from stars on the ZAMS (see for example Barrado 
y Navascues 1998) and therefore other age cali- 
brators have to be employed like the stellar mag- 
netic activity, the lithium abundance, the stellar 
rotation rate and the X-ray flux. In addition to 
these photometric signatures, the galactic space 
velocity is often employed for an age estimate and 
for defining young moving groups. Young stars 
are indeed characterized by low dispersion of the 
galactic space- velocity components (U, V, W). A 
few criteria exist in the literature defining pre- 
cise ranges for these components (Eggen 1996, Jef- 
fries 1995). These same criteria are employed by 
Montes et al. (2001). 

We avoided the selection of binary systems with 
separations smaller than ~ 2" because the proba- 
bility that a planet can find a stable orbit around 
one of the binary component is rather small. The 
binary systems included in our survey are: TWA 
5 (2"), TWA 8 (13"), TWA 9 (6"), GJ 799 (278) 
and BD -17° 6128 (272). 

3. NACO Observations 

The observations were performed with the AO 
system NACO, consisting of the infrared cam- 
era CONICA (Lcnzcn et al. 1998) assisted by the 
adaptive optics facilities NAOS (Rousset et al. 
2000) placed at the focus of the unit "Ycpun" of 
the ESO/VLT 8 m telescope. We observed 30 tar- 
gets during three observing runs: 16-19 February 
2003, 20-23 July 2003 and 7-9 January 2004. Two 
targets (GJ 182 and GJ 179) were observed in ser- 
vice mode. Table 3 summarizes the observational 
instrumental parameters. 

Due to the rather bright central star, our ob- 
servational strategy was the following. First, we 
observed each target in a narrow-band filter (Ta- 
ble 3) for a few minutes in order to obtain a non- 
saturated point spread function (PSF). We ob- 
tained a Strehl ratio (SR) in the range 0.2-0.55 
over the all runs. Second, we obtained deep and 
saturated (~ 20-30 min) K s or H broad-band ob- 
servations of each target (Table 3) in order to in- 
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crease the detectability on the PSF wings. We 
used NACO at a pixel scale of 13.25 mas yielding 
a FOV of 14" x 14". The optical wavefront sen- 
sor was locked on the target star yielding a typical 
resolution of ^50 mas in K s band. 

For the deep observations, we selected the ex- 
posure time (DIT in Table 3) of a single frame to 
obtain a number of counts (for the central star) 
larger than <~ 10 — 30 times the saturation limit. 
The deep images are saturated over a central di- 
ameter having a size of 0"l-0"2. For each target, a 
data set of dithered frames was obtained in order 
to correct for sky emission and the pixel to pixel 
variations in sensitivity. 

In a few cases (e.g. targets characterized by 
a smaller brightness) we used a broad-band filter 
with a neutral density filter (ND Short - see Ta- 
ble 3) instead of a narrow-band filter to obtain 
a non-saturated PSF. For HD 17925, due to its 
high brightness, we used a different observational 
strategy based on the absorption properties of a 
planet atmosphere. This target was observed with 
the intermediate-band filters (IB 2.09 and IB 2.24) 
centered out and in the methane (CH4) absorption 
band (~ 2.24 //m). We expect that planet features 
would appear in the continuum (IB 2.09) but not 
in the absorption band (IB 2.24). 

Due to the higher Strehl ratio (SR) reachable in 
the K a band compared to the H band, we gave pri- 
ority to the former band in order to optimize the 
probability to detect a low-mass companion. We 
chose the H band only for the oldest stars (100-200 
Myr), i.e. for planets having a theoretical colour 
(H-K) < — 1.3 mag. We note that, according to 
the COND models (Baraffe et al. 2003) (see Fig.l) 
and for planet masses of 3-10 Mj, the H band is 
comparable to the K band for an age minor than 
100 Myr. For targets with an age of 100-200 Myr 
the gain of the H band with respect to the K band 
is significant. 



4. Analysis 

4.1. Data processing 

The procedure used to reduce our data set in- 
cludes IRAF 4 and ECLIPSE 5 packages. Each 
frame was cleaned by a bad-pixel mask and di- 
vided by the flat field. In addition, the sky back- 
ground was eliminated in each frame. The fi- 
nal deep image (shift and add of all the frames) 
was obtained using a cross-correlation routine in 
ECLIPSE. We underline that the shift and add 
had to be done using a cross-correlation method 
since our PSFs are saturated in the central zone. 
Alternative methods, based on the gaussian fit to 
the peak of a PSF fail in such a case. 

In order to save observational time, we did not 
observe any reference standard star but we filtered 
the final deep images from the low spatial frequen- 
cies using two different methods. In this way the 
large-scale features of the saturated PSF were re- 
moved. 

In the first method, we convolved the deep im- 
age of the scientific target with a two-dimensional 
Gaussian function having a width roughly equiv- 
alent to the FWHM of the PSF of the star. Then 
we subtracted the result from the original deep 
image. Only features identified outside the cen- 
tral saturated region (~ 0"l-0"2) are reliable in all 
cases. We note that, if the FWHM of the Gaus- 
sian function is much smaller than the FWHM of 
the PSF of the star, the filtering is not efficient. 
This means that large-scale features of the satu- 
rated PSF are not sufficiently eliminated and the 
potential planet-like features close to the star are 
not visible. If, on the contrary, the FWHM of the 
Gaussian function is too large, the filtering elim- 
inates also the potential planet-like features close 
to the star i.e. the structures having small spatial 
scale. 

In the second method, we selected, for each 
point of the deep image, an annulus. Calculat- 
ing the average of these frames we obtain the esti- 
mate of the background of the deep image inside 
an annulus around the centre of our target. We 

4 IRAF is distributed by the National Optical Astronomy 

Observatories, http://iraf.noao.edu/ 
5 ECLIPSE is a data reduction package developed by ESO, 

N. Devillard, "The eclipse software", The messenger No 87 

- March 1997 
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finally subtracted the background from the orig- 
inal image. The annulus is characterized by an 
inner radius of half of the PSF FWHM and an 
outer radius ^3 larger than the internal radius. 
The central part of the PSF is excluded in the cal- 
culation of the background because it introduces 
an offset in the average estimation. Reducing the 
size of the inner radius, we increase the filtering of 
the low spatial frequencies. Increasing the outer 
radius the estimation of the background is affected 
by features belonging to regions far away from the 
central part of the PSF. We define an outer radius 
equal to 3 times the inner radius in order to esti- 
mate a background on a limited region around the 
central part of the PSF. We note that the final de- 
tection limit only weakly depends on the values of 
the inner and outer radius. Repeating the filtering 
with several combinations of inner and outer radii 
the final detection limit shows differences within 
0.2 mag. 

The filtering was applied in a region of ~ 6" x 6" 
centred on each target. Although the two methods 
have differences, they are based on the same prin- 
ciple aiming to eliminate the large-scale features 
of the saturated PSFs. In Section 4.3 the similar- 
ities/differences of the two methods are discussed 
with respect to the results that we obtained. 

4.2. Point-like features 

After data reduction and filtering of all the 
observations as described in Sec.4.1, a visual in- 
spection of all the deep images was performed and 
a few point-like structures were identified in some 
of the images: V343 Nor (K s band), SAO 252852 
(H band), GJ 179 (H band), BD +2° 1729 (H 
band), BD +1° 2447 (H band), GJ 183 (H band) 
and HD 221503 (H band). All the other targets 
did not reveal any point-like sources around them. 
We discuss here the nature of these objects try- 
ing to distinguish between artifacts (ghost images) 
and faint point sources at low galactic latitude. 

We have verified that all the targets observed 
in H band show three point-like features at the 
same position at about 1" from the star. Fig. 2 
shows deep images of two of these targets: BD 
+2° 1729 and BD +1° 2447. Three point-like fea- 
tures are visible (indicated by arrows and letters 
A, B and C) at the same relative angular separa- 
tion from the star. These features are probably 



'ghost' images produced by some imperfection of 
the broad-band filter or static aberration residuals 
that depend on the wavelength. This could justify 
the fact that they are visible only in H band. 

In only two deep images (V343 Nor and SAO 
252852 - Fig. 3) we observed clear point-like 
sources different from those created by the H fil- 
ter that we have just described. Fig. 3a shows 
the deep image of SAO 252852 (galactic latitude 
= -4.6°). An object at around 25-30 AU from 
the star is visible. The estimated magnitude of 
this object is ~17 mag in H band. Due to the 
low galactic latitude of SAO 252852, the proba- 
bility that this feature is a background star is not 
negligible. Fig. 3b shows the deep image of V343 
Nor (galactic latitude = —1.9°). Several sources 
were detected around this star. The estimated 
magnitude of all the sources is in the range 15-20 
mag. We note that a 5Mj planet orbiting around 
this 12 Myr old star would have ~ 16 mag in K 
band. Besides this, we observe that the galac- 
tic latitude of V343 Nor is relatively low (—1.9) 
and several among these point-like sources appear 
a few arcsec away from the central star. Thus 
the probability that these sources are background 
stars is quite high. We can not totally exclude 
that the faint sources around SAO 252852 and 
V343 Nor are gravitationally bound massive plan- 
ets. For this reason, a second epoch observation 
is necessary in order to clarify the nature of the 
detected point-like sources. Knowing that the 
proper motions (Hipparcos catalog, Perryman et 
al. 1997) of SAO 252852 are {p a = -202 mas/yr 
and ns= -270 mas/yr) and those of V343 Nor 
are (fJ, a = -52.87 ± 1.16 mas/yr and fj,s= -105.99 
± 0.98 mas/yr) we deduce that after one year a 
re-observation could throw light on the nature of 
these objects. A forthcoming run is planned to 
definitely verify if these objects are indeed back- 
ground stars. We remind the reader that an es- 
timation of the probability that a point source 
having a faint magnitude comparable to ours in a 
field of view of 2" is a background star was done 
for a galactic latitude of ~ +15° (Brandner et 
al. 2000). A probability of 70% was estimated 
in that case. Our point sources are placed at a 
lower galactic latitudes (—1.9° and —4.6°), and 
the point-sources are placed at around 2" from 
the central star (in the SAO 252852 case, Fig.3) 
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and at distances larger than 2" (in the V343 Nor 
case, Fig. 3). We conclude, that the probability 
that such sources are background stars is larger 
than 70%. 

We conclude, thus, that no very promising 
massive-planet candidate was identified in our sur- 
vey. 

4.3. Estimate of detection limits 

To estimate quantitatively the detection limit of 
our deep images as a function of the angular sep- 
aration from the central star, we used the follow- 
ing procedure. Along a radial direction we calcu- 
lated the standard deviation of the intensity over 
a box of dxd pixels (with typical values of d of 
4-6 pixels which corresponds to the FWHM of the 
non-saturated PSFs) and one pixel step. We then 
averaged over all radial directions (360°) and fi- 
nally we calculated the contrast AM at 5a with 
respect to the peak of the non-saturated PSF of 
the star. Since the central star is saturated in 
the broad-band images we had to use the fluxes 
measured in the non-saturated narrow-band filters 
(see Table 3) to calculate what peak flux our satu- 
rated broad-band images would have if they were 
not saturated. The non-saturated PSF had to be 
multiplied by a constant factor (A) to take into 
account the different transmission and width of 
the filters. If F 2 is the filter used to obtain the 
saturated PSF and F\ the one used to obtain a 
non-saturated PSF, A is given by: 

FWHM F2 Tr F2 DIT F2 
~ FWHMfi ' Tr F i ' DITfi 

where Tr is the filter transmission coefficient, 
FWHM is the full width half maximum of the 
filters and DIT is the integration time of each sin- 
gle frame (Table 3). 

In order to illustrate that our detection limit 
is reasonable, we scaled the non-saturated PSF of 
one of our targets (V343 Nor) to a flux value cor- 
responding to a 5Mj planet (representative of the 
typical planets we are interested in) and we placed 
this 'artificial planet' in the deep image at different 
angular separations from the star. Fig. 4 (a)-(b) 
shows the filtered deep images of the same target 
processed with the two filtering methods described 
above (Section 4.1). This "artificial planet" is 
placed at 1" and 077 from the central star. Fig. 4 



(a)-(b) shows that a separation of 0'.'7 is roughly 
the minimum angular separation at which a 5Mj 
planet can be detected. 

Figures 5-8 show the detection limits as a func- 
tion of angular separation for 28 objects of our sur- 
vey for which deep images have been obtained. We 
note that, in one case (GJ 207.1), a technical prob- 
lem did not allow to observe the non-saturated 
PSF and this prevented us to calculate the detec- 
tion limit. In the case of the binary GJ 799, we 
could calculate only the detection limit for GJ 799 
A. The separation of ~ 2'.'8 between the two com- 
ponents was slightly too large to include both stars 
in the same frame in all the dithered positions. 

For all targets, the detection limit is reliable 
starting for angular separations > O'.'l. For angu- 
lar separations < O'.'l the PSF of the deep images 
is usually saturated. Fig. 8 shows the detection 
limit of the two targets (HD 17925 and BD -17° 
6128) that we observed through the more narrow 
"methane band filters" (see Section 3). Since the 
optical quality of these intermediate-band filters 
is worse compared to the broad- and narrow-band 
filters, the final deep images appear dominated by 
structures characterized by high spatial frequen- 
cies. For this reason, the filtering is less efficient 
for these two deep images (the dashed, dotted and 
full thin lines are quite similar - Fig. 8). 

From Fig. 5, 6, 7 we can see that: 

• The detection limit of the filtered deep im- 
ages is very similar for the two different 
filtering methods. This indicates that the 
achieved gain from the applied 'filtering' is 
practically method independent. 

• The filtering of large spatial structures (low 
spatial frequencies) has a variable efficiency. 
The angular separation range over which the 
filtering shows its maximum efficiency is 0'.'5- 
1'.'5 from the central star. In this region the 
gain in contrast AM with respect to the non- 
filtered image ranges from 0.5 to 2 mag. 

• The detection limit is not the same for all 
the stars. This is not surprising because 
it depends on several parameters: the age, 
the magnitude of the central star, the wave- 
length used for the observations, the dis- 
tance of the star from the observer and also 
on the exposure time. In one case (the brown 



6 



dwarf TWA 5B), the detection limit is rather 
low because we did not saturate the PSF 
since the central star is relatively faint. 

To calculate the typical performances of NACO 
at subarcsecond distances, we averaged the con- 
trast AM obtained at 0'/5 and 1" for all the 
observed targets having a total integration time 
within the range 20-30 minutes. We find, in the 
K s band, a AM = 9 mag at 0'.'5 and a AM = 11.5 
mag at 1". We find, in H band, a AM = 9.2 mag 
at 0'.'5 and a AM = 11.7 mag at 1". We recall that 
0'.'5 corresponds, for our sample, to projected sep- 
arations of 2-35 AU, 1" to projected separations 
of 4-70 AU. 

Luhman & Jayawardhana (2002), observed 25 
targets with NICMOS (Keck II) and found a AM 
~ 10-10.5 mag at 1" in the H band with an in- 
tegration time of ~ 18 minutes. This illustrates 
nicely that NACO is a well performing instrument 
for high-contrast imaging observations. 

5. Discussion 

5.1. Upper limit of the projected separa- 
tion between the star and a potential 
planet 

The detection limit estimates also permit us to 
calculate an upper limit of the projected separa- 
tion between a star and a potential planet, as fol- 
lows. We consider a typical 5Mj planet, we calcu- 
late the theoretical value (retrieved by the mod- 
els) of the star/planet contrast, we then place this 
value on the ordinates of the detection limit curve 
(see for example Fig. 5) and finally we retrieve (on 
the abscissae) the distance beyond which we can 
affirm that such a planet should not exist (because 
it was not detected). Table 4 shows, for each tar- 
get, such upper limits for the star/exo-planet dis- 
tance. Table 5 shows the frequency at which we 
would be able to detect a 5Mj planet at different 
distances from the parent star. By analysing these 
results from a statistical point of view, we can cal- 
culate (Fig. 9) the cumulative distribution of the 
upper limit of the distance star/exo-planet for our 
sample of targets. The thin line represents the 
cumulative distribution obtained for 5Mj planets. 
We find that (Fig. 9), with respect to our sample of 
targets, in 50% of cases (median value), there are 
not planets at distances larger than 14 AU and, 



in 100% of cases, there are no planets at distances 
larger than 65 AU. We underline that we are con- 
sidering a 5Mj planet just as a reference, but the 
same calculation could be done for more and/or 
less massive planets in the range 3-10 Mj. In the 
case of more massive planets, the cumulative dis- 
tribution should have a median value (50%) asso- 
ciated with an angular separation smaller than 14 
AU. In the case of less massive planets, the cumu- 
lative distribution should have a median value (50 
%) associated with an angular separation larger 
than 14 AU. To better appreciate the sensitivity 
with respect to the mass of the planets we also 
calculated the cumulative distribution for 10Mj 
planets (Fig.9 - bold line). We find that, in 50% 
of the cases (median value) there are no planets at 
distances larger than 8.5 AU and in 100% of the 
cases there are not planets at distances larger than 
36 AU. Table 4 and Table 5 show the upper limit 
of the projected star/planet separation and statis- 
tical results also for 10Mj planets. The numbers 
given in these tables illustrate that it is relatively 
unlikely that a massive planet exists beyond a cer- 
tain distance from the star. For example the fact 
that no 5Mj planet was found for 18 of our tar- 
gets beyond 20 AU (see Table4) means that the 
likelood for the existance of such planets beyond 
this distance is - (1/18) % i.e. - 5.6%. This 
last estimation is useful for the reader to have an 
idea of how much the size of the sample affects 
the statistical results. We recall that we are using 
the term probability in a less conventional way. A 
precise (statistically speaking) definition of proba- 
bility should need some positive detection that we 
obviously do not have. 

Which are the implications of these results with 
respect to models such as those cited in Section 
1 (gravitational scattering model and outward mi- 
gration) , explaining the transfer of planets to large 
distances (some tens of AU) from the central star 
? Our results indicate that the frequency of planet 
'jumps' or outward migration is quite low. More 
precisely, it seems difficult for 10Mj planets to find 
a stable orbit at distances larger than 16 AU, and 
that 5Mj planets at distances larger than 26 AU, 
at least in the range of early ages that we consid- 
ered. 

We note that the statistical analysis presented 
in Fig.9 related to the 10Mj planets is particularly 
interesting because, due to a steeper cumulative 
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distribution concentrated at small projected sepa- 
rations from the central star, it permits us to bet- 
ter constrain the probability of location of planets. 
We also recall that such massive planets are rarely 
detected at distances smaller than 6 AU using ra- 
dial velocity measurements. As described in the 
introduction, the radial velocity technique is sen- 
sitive only to planets orbiting at closer distances 
(< 6 AU) from the central star. It is interesting to 
analyse what the direct imaging technique reveals 
about lOMj planets at distances larger than 6 AU. 
Looking at our results, we can state that, for most 
of the targets (~ 90% of cases), there are no lOMj 
planets at distances larger than ~ 16 AU. We find 
that the probability to find planets equal or more 
massive than lOMj at large distances is low as well 
as the probability to find them at close distances 
(< 6 AU) as appears from radial velocity surveys 
(Marcy et al. 2003). Massive planets are expected 
to form beyond the so called 'snow line ' distance 
(~ 5 AU for the solar-like stars 6 ) (Lissauer, 1993) 
and from this distance they might migrate inward. 
Our study shows that, at least for late spectral 
type (K and M) stars, if such massive (~ lOMj) 
planets form beyond the snow line, we can envis- 
age the following possibilities: (1) massive planets 
form and can be found in the narrow range of dis- 
tances 'snow line ' - 16 AU from the central star, 
(2) at distances larger than 16 AU only planets 
having masses smaller than lOMj can form and 
can be found, (3) if any lOMj planet would have 
formed beyond the 'snow line ', the formation and 
inward migration necessarily happen on a time 
scale < lOMyr. To test the hypotheses (1) and (2) 
one should improve the contrast at small angular 
separations. To test the hypothesis (3) one should 
select a sample of stars younger than ours. Un- 
fortunately, stars having this property are placed 
preferably at distances > 100 pc from the observer. 
Even future ground-based planet finders will be 
able to detect planets only at larger separations (a 
few tens of AU) from the central star if the targets 
are placed at such distances (> 100 pc). For this 
reason, it would be probably more efficient and 
appropriate to observe them from space. In this 
way, one could by-pass the limitations imposed by 

'We underline that the 'snow line' is equal to 5 AU for 
typical solar-like star. This distance is calculated for simple 
thermal equilibrium and is estimated equal to 1-2 AU for 
late type stars (Black 1980) 



the atmosphere on the attainable angular resolu- 
tion and to probe the same projected separations 
(as those analysed in this paper) from these dis- 
tant stars. The resolution reachable in space by 
the Hubble Space Telescope (HST) is not good 
enough to probe such small physical distances. In 
a previous study (Brandncr et al. 2000), a set of 
pre-main sequence T Tauri stars in the Chamaleon 
T (1-5 Myr at 150 pc) and Scorpius-Centaurus OB 
(5-15 Myr at 150 pc) associations were observed. 
Results showed that at 30 AU it would be possi- 
ble to detect only more massive planets (~ 20Mj) 
than those we would detect with NACO (this pa- 
per) . We should thus wait for a new generation of 
space telescope to probe such close distances at so 
early ages. 

Other interesting issues can be retrieved by our 
analysis concerning the core-accretion model. Our 
results provide interesting information concerning 
mechanisms to explain planet formation. We are 
considering here the core-accretion model origi- 
nally proposed by Pollack (1996) since it provides 
an explanation for the formation of all gaseous as 
well as solid planets 7 . It claims that a growing 
solid core, after reaching a critical mass, accretes a 
massive atmosphere. According to our results, this 
model could work, assuming that massive plan- 
ets (> lOMj) form in a narrow 5-15 AU distance 
from the central star in a time scale of ^10 6 yr 
(typical time scale of disk lifetime). One of the 
critical aspects of the core-accretion model has al- 
ways been that it can not explain the presence 
of planets observed with the radial velocity tech- 
nique at close distances from the central star. The 
time necessary for the formation plus the migra- 
tion should be longer than the typical disk lifetime 
(-10 6 yr, Haisch et al. 2001). The modified ver- 
sion of the Pollack model, recently proposed by 
Alibert et al. (2004), seems to solve this prob- 
lem. Using a detailed description of the compo- 
sition (gas plus solid) of the disk this model re- 
duces the time necessary for a planet to form and 
migrate by a factor of 10. With such a modifica- 
tion, the process of formation and migration might 
need a time shorter than the typical disk lifetime. 
Following the predictions of Alibert 's model, it 
would be possible that we did not detect mas- 



We precise that the core-accretion model was proposed 
for solar-like star and our sample is done by young stars 
havinging later spectral type than solar-like stars. 
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sive planets simply because they formed at ages 
(< lOMyr) earlier than those of our selected sam- 
ple. This means that our results are consistent 
with the core-accretion model in the Alibert for- 
mulation. This model seems thus promising due 
to the fact that it matches with both observational 
issues coming from radial velocity measurements 
and also from direct imaging observations. 

We finally mention that our results indicate 
that a sort of 'exo-planet desert' for massive (~ 
lOMj) planets is observed at distances larger than 
15 AU from the stars of our survey. This desert 
seems to extend to quite large distances from the 
central star. A recent survey (McCarthy & Zuck- 
erman, 2004) carried out with the near-infrared 
camera NIRC and a coronograph at Keck did not 
detect any 5-12 Mj exo-planet in the range 75-300 
AU around ~ 42 stars having a mean age of 300 
Myr. 

5.2. Level of sensitivity: comparison with 
previous surveys 

We would like to mention the relevant improve- 
ment in terms of sensitivity reached by our study 
(using NACO at the VLT) with respect to previ- 
ous published results. A recent analysis of massive 
planet detection (Neuhauser et al. 2003), carried 
out with the MPE Speckle Camera at the ESO's 
3.5 m NTT telescope, put an upper limit of 60 
AU on the distance between the star and 10Mj 
exo-planets in the case of targets belonging to /3 
Pictoris (~12 Myr). Results from our study show 
that we can exclude any 10Mj planets at distances 
larger than 15 AU (a physical distance 4 times 
smaller) for targets belonging to the same mov- 
ing group. Our results show a better sensititity 
also with respect to previous observations made 
with HST (Neuhauser et al. 2002). Estimates pro- 
vided by this study give an upper limit distance 
between the star and 5Mj planets of ~ 50 AU. 
Stars observed by Neuhauser et al. (2002) belong 
to a ~2Myr old systems (well younger than our 
sample - 10-200 Myr). Models predict a mean dif- 
ference in H and K band of the order of at least 
3 mag between the age of stars in Neuhauser et 
al. study and ours. We thus deduce that the same 
instrument (HST) would have an upper limit for 
the star/planet distance quite a bit larger than 50 
AU for 5Mj planets in the 10-200 Myr age range. 



5.3. Forthcoming observational strategy 

How can we improve our results to better con- 
strain planet-formation models ? It is evident that 
an increase of the achievable contrast at smaller 
angular separations shifts the median value of 
the cumulative distributions to smaller angular 
separations (Fig. 9). This means that the upper 
limit of the projected star/exo-planet separation 
decreases. However, we note that, using the ob- 
servational strategy employed in our survey, we 
can hardly improve the sensitivity at subarcsec- 
ond distances using a 8 meter class telescope as 
illustrated in Fig. 10. This figure shows the de- 
tection limit obtained for one of our targets for 
an integration time of 5, 10 and 22 minutes. The 
deep image was filtered from low spatial frequen- 
cies before we calculated the detection limit. We 
observe that, beyond 1", the reachable contrast 
increases with increasing the integration time (~ 
1 mag in 15 min). At distances smaller than 1", 
on the contrary, a longer integration time does 
not increase the contrast because we are strongly 
limited by the speckle noise. 

The differential imaging technique proposed by 
Racine et al. (1999) is a rather efficient method to 
reduce the speckle noise at small angular separa- 
tions from the central star. Images taken simulta- 
neously at slightly different wavelengths are sub- 
tracted from each other. Since the speckle noise 
of both images is nearly the same 8 , the speckle 
noise is eliminated in the resulting image provid- 
ing one can sufficiently limit static aberrations. 
The differential technique finds a useful applica- 
tion in detecting cool methane-rich faint objects. 
Rosenthal et al. (1996) first proposed observa- 
tions of the same star done simultaneously in two 
different wavelengths which are located "on" and 
"off" the methane (CH4) absorption band. If one 
subtracts the image taken in the absorption band 
from the one taken outside the absorption band, 
planet/brown dwarf features appear in the con- 
tinuum. The absorption at 1.62 fim, in the CH 4 
band, is a clear signature of any object having 
temperature T < 1300 K (giant exo-planets and 
T-type brown dwarfs). 

We precise that speckles of images taken at slightly different 
wavelengths match if frames are rescaled opportunely by a 
A2/A1 factor. 
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We can thus increase the sensitivity (contrast) 
at subarcsecond distances using the recently im- 
plemented NACO/SDI (Simultaneous Differential 
Imager) instrument. This instrument employs the 
differential technique in H band. Preliminary re- 
sults (Close et al. 2004, Lenzen et al. 2004) indi- 
cate a AM ^11 mag at Of! 5 from the central star. 
Considering our averaged estimate of AM = 9.2 
mag in the H band (see Section 4.3), we should 
gain ~ 2 mag in contrast using NACO/SDI at Of! 5 
(see Fig. 10). 

A further and fundamental improvement in 
contrast should be attained by future ground- 
based "planet finders" , i.e. instruments supported 
by AO systems providing a better Strehl Ratio. 
The contrast reachable by such an instrument at 
0'/5 should be AM ~ 17.5 mag in J-H bands (see 
Fig. 10). We refer in particular to the top level 
requirements of CHEOPS, a second generation in- 
strument for the VLT conceived for planet detec- 
tion (Gratton et al. 2004). We think that the 
principal goals of these instruments (NACO/SDI 
and ground-based planet finders 9 ) in the context 
of the detection of massive and young exo-planets 
and the understanding of mechanisms responsible 
for exo-planet formation are the following: (A) 
They will permit us to provide answers to the hy- 
pothesis (1) and (2) (see Section 5.1). (B) They 
will be rather appropriate to do studies similar to 
the one proposed in this paper but applied to a 
sample of targets characterized by an earlier spec- 
tral type (i.e. solar-like type) in the same range 
of ages. This is due to the better contrast reach- 
able at subarcsecond distances (see Fig. 10). (C) 
They will be appropriate instruments to detect 
even fainter exo-planets (i.e. in the 1-5 Mj range) 
at subarcsecond angular separations from the cen- 
tral star. 

6. Conclusion 

In this paper we have presented results of a sur- 
vey carried out with NACO at the VLT aiming 

9 We underline that, in the context of this paper, we are 
talking about the category of planets in which the intrin- 
sic luminosity is quantitatively dominant over the reflected 
one. Planet finders will also be able to detect colder and 
older exo-planets for which the reflected light is dominant 
with respect to the intrinsic one. This other category of 
planets orbit in general at closer distances from the central 
star. 



to detect massive exo-planets (3-10 Mj) orbiting 
around 30 young (10-200 Myr) and nearby (d < 
77 pc with 70 % of them in the first 32 pc) stars. 

No promising candidate was identified. All faint 
sources found in the vicinity of the central star are 
presumably background stars. A detailed estimate 
of the detection limit in 28 of the observed targets 
was presented. This permitted us to calculate the 
typical contrast (AM) between a planet and a cen- 
tral star at C/5 and 1" provided by NACO for a 
typical integration time of 20-30 min (see Section 
4.3), to present a statistical analysis of the upper 
limit for the star/planet separation and also to 
discuss our results in the framework of the mecha- 
nisms that are supposed to explain the formation 
of exo-planets. 

Our most relevant result is that, in our sample 
of targets, in 100% of cases, no 10Mj (and more 
massive) exo-planets was detected at distances 
larger than 36 AU and no 5Mj (and more mas- 
sive) exo-planets was detected at distances larger 
than 65 AU. If one considers the median values of 
the cumulative distributions we have that, in 50% 
of cases, no 10Mj (and more massive) exo-planets 
should exist at distances larger than 8.5 AU and 
no 5Mj (and more massive) exo-planets should ex- 
ist at distances larger than 14 AU. We emphasize 
that, these statistical estimates were done for the 
projected star /exo-planet distance. We are not 
considering here the possibility that planets might 
not be visible in part of their orbit. Our study in- 
dicates that massive planets (mass > 5Mj) are 
rare at distances larger than 6 AU as well as at 
distances smaller than 6 AU as indicated by ra- 
dial velocity estimates (see Section 5.1). 
We underline that for the category of massive, 
warm and self-luminous planets that we looked for 
in our survey, it is not possible to calculate a re- 
liable probability distribution for detecting plan- 
ets. The reason is that, to do this, one should 
need to know the planet frequency distribution 
'a posteriori'. This is obviously impossible to 
do if no-planets have already been detected. For 
this category of planets, one can only define some 
merit functions (depending on several parameters 
such as the age, the distance from the observer, 
the magnitude of the central star, the mass of 
the planet) to define some criteria to select tar- 
gets for which the detection should be particularly 
favourable. For this reason, we think that it is fun- 
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damental to report results (even of non-detections) 
for the category of exo-planets that we are search- 
ing for and to try to get as homogeneous analysis 
as possible of the results in order to improve the 
statistics of the upper limit for star/planet dis- 
tance. 

7. Appendix 

GJ 803 (Au Mic): 

Recent estimates of emission at optical wave- 
lengths (Liu ct al. 2004, Kalas et al. 2004) of a 
circumstellar disk surrounding GJ 803 (AU Mic) 
were published. The disk is detected between ~ 
50 and 210 AU corresponding to around 5" - 21" 
(d ~ 10 pc). The inner region was occulted by 
the presence of a coronograph so no information is 
provided on the disk structure at distances smaller 
than 5". The field of view of our observations is 
much more smaller (~ 2" in radius around the 
central star) than that used for the disk detection. 
Within ~ 50 AU, the typical time scale of dust 
is shorter than the stellar age so it was suggested 
(Kalas et al. 2004) that grains observed here must 
be continually replenished by the collisional ero- 
sion of much larger objects (like comets and aster- 
oids). In the inner region of the disk, the presence 
of planetisimals might originate planets by accre- 
tion. Kalas et al. (2004) suggest that in the 2-30 
AU region, the presence of gas giant planets should 
be possible. Our results indicate that no planets 
equal or more massive than 5Mj are detected at 
distances larger than 4 AU and no planets equal 
or more massive than 10M,/ are detected at dis- 
tances larger than 2 AU. Besides this, we can not 
exclude that some planets with mass smaller than 
5Mj exist at such close angular separations. We 
have to wait for observations at higher contrast to 
have a more precise answer to this question. Also 
it was suggested (Liu et al. 2004) that the spa- 
tially localized enhancements and deficits found 
in the GJ 803 disk at distances larger than 15 AU 
from the star might be signposts of ongoing planet 
formation. Our observations can not infer any in- 
formation at this large distance from the star. 

We thank the ESO staff for their support during 
the observations. We also acknowledge I. Baraffc 
and A. Burrows for providing informations on 
their models and W. Brandner and A. Raga for 
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Fig. 1. — Theoretical absolute magnitudes in the 
K and H bands vs. planet mass (in Jupiter mass 
units) according to the COND models of Baraffe 
et al. (2003). The light lines (green in on-line 
version) refer to the H band and the dark lines 
(black in on-line version) to the K band. The full 
thin line represents an age of 10 Myr, the dotted 
thin line an age of 30 Myr, the dot-dashed bold 
line an age of 100 Myr. 
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Fig. 2.— Two deep images obtained through the H band filter (BD +2° 1729 on the left, BD +1° 2447 
on the right). Both images have been processed using a high pass spatial frequency filter. Three point-like 
sources (A, B and C) are visible at the same position, at around 1" from the central star in both images 
and in all other deep H band images. This clearly indicates that these point-like sources are ghost images 
produced by the H band filter. 
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Fig. 3. — Left: a point source (marked with a circle) is visible at around 30 AU from the star SAO 252852 
(galactic latitude -4.6°). The deep image was processed using a high pass spatial frequency filter. A, B and 
C indicate the ghost images due to the H band filter (see Fig. 2). Right: several point sources (marked with 
circles) are detected at several tens of AU from the star V343 Nor (galactic latitude -1.9°). In this picture 
the deep image is not processed with a high pass spatial frequency filter because all point sources appear at 
distances larger than 2" from the central star. The encircled point sources are probably background stars 
since both stars are at low galactic latitudes. 
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Fig. 4. — (a) Deep image of V343 Nor (2'.'7 x 2'.'7) obtained in the Ks band (~ 22 min) and processed with 
a high pass spatial frequency filter. The PSFs of two 5Mj planets ('artificial planets') were obtained by 
scaling a non-saturated PSF of the central star. Two artificial planets were placed at 0'.'7 and 1" from the 
central star, (b) The same as (a) but using a different filtering method, (c) Contrast expressed in AM (at 
5cr) vs. angular separation derived from the filtered images shown in (a) and (b). The dashed line represents 
the detection limit of the non-filtered PSF of the star, the dotted and thin full lines represent the detection 
limit of the PSF of the star filtering out low spatial frequencies by using the two different methods shown 
in (a) and (b) (see text). The x-axis projection of the black dot represents the minimum angular separation 
at which a 5Mj planet should be detectable. The 077 value corresponds to the visible planet in the (a) and 
(b) pictures. The grey scale is arbitrary, the contrast in the picture (b) is better than in picture (a). 
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Fig. 5. — Planet detection limit (or contrast) expressed in AM (at 5er) vs. angular separation for 9 observed 
stars. The dashed line refers to the non-filtered images, while the full thin line and dotted line refer to the 
images processed by a high spatial frequency filter using two different methods (see text). ET represents the 
total exposure time. The exposure time of the single frame (DIT) is given in Table 3. No particular criterion 
is used for the order of the targets in this figure as well as in Fig. 6, Fig. 7 and Fig. 8. 



18 



Fig. 6. — As Fig.5 for further 9 stars. The bumps visible in GJ 799A and BD -17°6128 corresponds to the 
companion of the binary system. 
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Fig. 7. — As Fig. 5 for further 9 stars. 



Fig. 8. — As Fig. 5 for further 2 targets. 
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Fig. 9. — Cumulative distribution of the upper 
limit for the projected star/planet separation cal- 
culated with our sample of targets. Thin line: 
5Mj, bold line: lOMj. This distribution implies 
that in 50% of cases (median value) we have not 
found any 5Mj planets at distances larger than 14 
AU and no lOMj planets at distances larger than 
8.5 AU. 
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Fig. 10. — Detection limit vs. angular separation 
of one target of the survey obtained with integra- 
tion times of 5, 10 and 22 min. The two horizon- 
tal bars placed at 0'.'5 mark the expected contrast 
reachable with SDI/NACO and CHEOPS (VLT 
planet finder - see text). 
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Table 1 

poperties of the central star^ 



Targets 


a (2000) 


5 (2000) 


ST 


D 


V 


Group Ass. 


Age 


Gal. Lat. 


Refer. f b > 




(h m s) 


(o ' ") 




pc 


(mag) 




(Myr) 


(°) 




Hip 2729 


00 34 51.2 


-61 54 58 


K5V 


49.5 


9.56 


Tuc. Ass. 


35 


-55 


(5), (6), (7) 


HD 17925 


02 52 32.1 


-12 46 11 


Kl 


10.4 


6.0 


local ass. 


50 


-58 


(8), (12) 




04 52 05.7 


06 28 35 














( 2 ) 


GJ 182 


04 59 34.8 


01 47 01 


M0. 5 


26.7 


10.1 




35 


-24 


(2) . (10) 


Hip 23309 


05 00 47. 1 


-57 15 26 


K7V 


26.3 


10.09 


Beta Pic 


12 


-37 


csi no") n 31 


GJ 207.1 


05 33 44.8 


01 56 43 


M2.5 


16.8 


11.5 




50 


_2 


(1),(2) 


AO Men 


06 18 28.2 


-72 02 42 


K5V 


38.5 


10.09 


Beta Pic. 


12 


-28 


(31 (10) (131 


BD +2° 1729 


07 39 23.0 


02 01 01 


K7 


9.8 


9.82 


local ass. 


100 


12 


(2) 


LQ Hya 


09 32 25.6 


-11 11 05 


K0V 


18.3 


7.80 




50 


28 


(1),(8) 


TWA 6A 


10 18 28.8 


-31 50 02 


K7 


77.0 


12.0 


TW Hya. 


10 


21 


(4), (5) 


BD 1° 2447 


10 28 55.5 


00 50 28 


M2 


7.2 


9.63 


local ass. 


100 


47 


(1),(2) 


TWA 5B (2 ff ) 


11 31 55.4 


-34 36 27 


M8.5 


50.0 




TW Hya. 


10 


25 


(4),(5) 


TWA SA M 1 




OR K 1 K K 
-iD Ol DO 






(R 11 2) 


± vv nya. 








TWA SB 


11 32 41.5 


-26 51 55 


M5 


21.0 


(R=13.85) 


TW Hya. 


10 


33 


(4), (5) 


TWA 9A( C ) (9") 


11 48 24.2 


-37 28 49 


Ml 


64.0< d ) 


11.3 


TW Hya. 


10 


24 


(4), (5) 


TWA 9B 


11 48 24.2 


-37 28 49 


Ml 


64.0< d ) 


(R.= 12.98) 


TW Hya. 


10 


24 


(4), (5) 


SAO 252852 


14 42 28.1 


-64 58 43 


K5V 


16.4 


8.50 




100 


-4 


(2), (9) 


V343 Nor 


15 38 57.6 


-57 42 27 


K0V 


39.8 


8.14 


Beta Pic. 


12 


-2 


(3), (10), (13) 


V2306 Oph 


16 30 18.1 


-12 39 45 


M3.5 


4.4 


10.12 




100 


24 


(2) 


HD155555 AB* e ' 


17 17 25.5 


-66 57 02 


G5IV 


31.4 


7.20 


Beta Pic. 


-16 


24 


(3), (10), (13) 


HD 155555 C 


17 17 25.5 


-66 57 00 


M4.5 


31.4 


12.71 


Beta Pic. 


12 


-16 


(3), (10), (13) 


CD -64° 1208 


18 45 36.9 


-64 51 48 


M0 


29.2 


9.20 


Beta Pic. 


12 


-24 


(3), (10), (13) 


PZ Tel 


18 53 05.9 


-50 10 50 


K0 


49.6 


8.42 


Beta Pic. /Tuc. Ass. 


12 


-21 


(3) 


GJ 799A (2 ;/ 8) 


20 41 51.1 


-32 26 07 


M4.5 


10.2 


11.0 


Beta Pic. 


12 


-36 


(3), (10), (13) 


GJ 799B 


20 41 51.1 


-32 26 09 


M4.5 


10.2 


11.0 


Beta Pic. 


12 


-36 


(3), (10), (13) 


GJ 803 (Au Mic) 


20 45 09.3 


-31 20 24 


M0 


9.9 


8.61 


Beta Pic. 


12 


-37 


(3), (10), (13) 


BD -17° 6128 (2"2) 


20 56 02.7 


-17 10 54 


K7 


47.7 


10.60 


Beta Pic. 


12 


-35 


(3), (10), (13) 


GJ 813 


20 57 25.4 


22 21 45 


M2 


13.6 


12.0 




200 


-15 


(2) 


GJ 890 


23 08 19.5 


-15 24 35 


M0 


21.8 


10.84 




100 


-63 


(1),(2),(H) 


HD 221503 


23 32 49.4 


-16 50 44 


K5 


13.9 


8.60 


local ass. 


> 100-200 


-69 


(1),(2),(9) 



( a )egend, from the left to the right: name of the target, right ascension, declination, spectral type, distance of the target, magnitude in V band, name 
of the association or group to which the target belongs, age, galactic latitude, references 

(^Ref. related to the distance and age of the central star: (1) Barrado y Navascucs priv. comun.; (2) Montcs ct al. 2001; (3) Zuckcrman ct al. 2001a; 
(4) Webb ct al. 1999; (5) Song ct al. 2003; (6) Zuckerman &; Webb 2000; (7) Zuckerman et al. 2001b; (8) Wichmann et al. (2003); (9) Povcda ct al. 
(1994); (10) Barrado y Navascucs et al. (1999); (11) Barrado y Navascucs 1998; (12) Favata ct al. (1998); (13) Barrado y Navascucs 2001 

( c ) Based on their location in IR color-magnitude or HR diagrams, both TWA9A and TWA9B might be older than other members of the TWA association 
whose parallaxes have been measured by Hipparcos (Barrado y Navascues, Mohanty &c Jayawardhana (2004)) 

( d ^A mean value of trigonometric and photometric distance was used (sec (5)) 

( e ) Spectroscopic binary (Macintosh ct al. 2001) 



Table 2 

Infrared magnitude for the central star and the possible planetary companion^ 



Targets 


(b) 


™>i!h 




™(5Mj,H) 


A(M K ) 


A(M H ) 


(H-K) 5Mj 




(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


Hip 2729 


5.96 


6.13 


18.60 


18.75 


12.64 


12.62 


-0.02 


HD 17925 


3.50 


3.60 


16.28 


15.99 


12.78 


12.39 


-0.4 


GJ 179 


7.77 


8.01 


20.24 


18.31 


12.47 


10.30 


-2.17 


GJ 182 


6.49 


6.66 


17.25 


17.40 


10.76 


10.74 


-0.02 


Hip 23309 


6.37 


6.54 


15.19 


15.59 


8.82 


9.05 


0.23 


GJ 207.1 


7.20 


7.41 


17.33 


17.04 


10.13 


9.63 


-0.29 


AO Men 


7.80 


7.97 


15.98 


16.38 


8.18 


8.41 


0.23 


BD +2° 1729 


6.61 


6.76 


18.11 


16.98 


11.50 


10.22 


-1.28 


LQ Hya 


5.97 


6.07 


17.51 


17.23 


11.54 


11.16 


-0.38 


TWA 6A (c) 


7.97 


8.17 


17.35 


17.65 


9.38 


9.48 


0.15 


BD 1° 2447 


5.58 


5.78 


17.44 


16.31 


11.86 


10.53 


-1.33 


TWA 5B (2") 


11.5 


11.95 


16.41 


16.71 


4.91 


4.76 


-0.15 


TWA 8 A 


7 44 

f .41 


7 79 


1 4 ^ 


1 4 

14.00 


7 no 


7 1 1 


0.14 


TWA 8B 


9.01 


9.36 


14.53 


14.83 


5.47 


5.52 


0.00 


TWA 9 A (9") 


7.68 


7.95 


16.95 


17.25 


9.27 


9.30 


0.1 


TWA 9B 


9.14 


9.34 


16.95 


17.25 


7.81 


7.91 


0.1 


SAO 252852 


5.68 


5.85 


19.22 


18.09 


13.54 


12.24 


-1.3 


V343 Nor 


6.31 


6.41 


16.05 


16.46 


9.74 


10.05 


0.31 


V2306 Oph 


5.93 


6.17 


16.30 


15.17 


10.37 


9.00 


-1.37 


HD 155555 AB 


5.62 


5.7 


15.54 


15.94 


9.92 


10.24 


0.32 


HD 155555 C 


7.25 


7.51 


15.54 


15.94 


8.29 


8.43 


0.14 


CD -64° 1208 


5.73 


5.88 


15.38 


15.78 


9.65 


9.90 


0.25 


PZ Tel 


6.69 


6.79 


16.54 


16.94 


9.85 


10.15 


0.3 


GJ 799A (27 8) 


5.54 


5.8 


13.10 


13.50 


7.56 


7.70 


0.14 


GJ 799B 


5.54 


5.8 


13.10 


13.50 


7.56 


7.70 


0.14 


GJ 803 (Au Mic) 


5.00 


5.15 


13.03 


13.44 


8.03 


8.29 


0.26 


BD -17° 6128 (2"2) 


7.13 


7.28 


16.45 


16.85 


9.32 


9.57 


0.25 


GJ 813 


7.95 


8.15 


20.50 


18.57 


12.55 


10.42 


-2.13 


GJ 890 


7.37 


7.52 


19.84 


18.71 


12.47 


11.19 


-1.13 


HD 221503 


5.78 


5.95 


> 20.55 


> 18.62 


>14.77 


>12.67 


> -2.10 



















legend, from the left to the right: name of the target, magnitude of the central star in K band {m s ,K ) 
and in H band (m s> jf), magnitude of a potential planet (retrieved from COND models) having a 5M,/ 
mass in K band (m Pj /f) and in H band (ra Pj ir), contrast star/planet in K and in H band, color (H-K). 

^'Calculated from V and spectral type 

^The magnitudes in H and K of all the targets belonging to TWA Hydrae refer to Webb et al. (1999). 
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Table 3 

Observational and instrumental parameters'") 



Targets 


Narrow-band 


Broad-band 


DIT. - Broad-band 

(sec) 


Exp. Time - Broad-band 
(min) 


Date 
(dd/mm/yyyy) 


Hip 2729 


NB 2.12 


K s 


1.30 


~ 28 


23/07/2003 


HD 17925 


IB 2.09 


IB 2.24 


1.70 


~ 23 


22/07/2003 


GJ 179 


NB 1.64 


H 


1.74 


~ 36 


26/12/2003 


GJ 182 


NB 2.12 


K s 


1.30 


~ 35 


26/12/2003 


Hip 23309 


NB 2.12 


K s 


2.06 


~ 20 


23/07/2003 


GJ 207.1 


K,/ND Short 


K s 


9.00 


~ 30 


08/01/2004 


AO Men 


NB 2.12 


K s 


2.50 


~ 20 


17/02/2003 


BD +2° 1729 


NB 1.64 


H 


1.00 


~ 17 


19/02/2003 


LQ Hya 


NB 2.12 


K s 


0.80 


~ 23 


09/01/2004 


TWA 6A 


NB 2.12 


K s 


5.00 


~ 14 


09/01/2004 


BD 1° 2447 


H/ND Short 


H 


0.36 


~ 21 


09/01/2004 


TWA 5B (2") 


NB 2.12 


K s 


5.00 


~ 31 


17/02/2003 


TWA 8A (13") 


NB 2.12 


K s 


4.20 


~ 19 


17/02/2003 


TWA 8B 


K,/ND Short 


K s 


19.00 


~ 22 


17/02/2003 


TWA 9A (9") 


K,/ND Short 


K s 


21.00 


~ 25 


18/02/2003 


TWA 9B 


if„/ND Short 


K s 


21.00 


~ 25 


18/02/2003 


SAO 252852 


NB 1.64 


H 


0.40 


~ 24 


21/07/2003 


V343 Nor 


NB 2.12 


K s 


1.90 


~ 22 


19/02/2003 


V2306 Oph 


NB 2.12 


K s 


1.30 


~ 25 


19/02/2003 


HD155555 AB 


NB 2.12 


K s 


1.20 


~ 8 


18/02/2003 


HD 155555 C 


NB 2.12 


K s 


5.40 


~ 38 


20/07/2003 


CD -64° 1208 


NB 2.12 


K s 


1.10 


~ 24 


22/07/2003 


PZ Tel 


NB 2.12 


K s 


1.56 


~ 30 


23/07/2003 


GJ 799A (2"8) 


NB 2.12 


K s 


0.58 


~ 23 


21/07/2003 


GJ 799B 


NB 2.12 


K s 


0.58 


~ 23 


21/07/2003 


GJ 803 (Au Mic) 
BD -17° 6128 (2"2) 


Ks/ND Short 


K s 


0.36 


~ 30 


21/07/2003 


NB 2.12 


K s 


2.50 


~ 27 


21/07/2003 


GJ 813 


NB 1.64 


H 


0.70 


~ 28 


23/07/2003 


GJ 890 


NB 2.12 


K s 


3.00 


~ 24 


21-22/07/2003 


HD 221503 


NB 1.64 


H 


0.70 


~ 22 


22/07/2003 



^egend, from the left to the right: name of the target, narrow-band filter, broadband-band filter, exposure time of the single 
frame, total exposure time, observations date. 
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Table 4 

Upper limit of the star/planet projected separation in the cases of a 5M, 7 and a IOMj 

PLANET. 



Targets 


Age 


D 


proj. sep.-AU 


proj. sep.-AU 




(Myr) 


(pc) 


(5M. 7 ) 


(IOMj) 


Hip 2729 


35 


49.5 


64 


36 


GJ 179 


200 


12.1 


9 


5 


GJ 182 


35 


26.7 


19 


11 


Hip 23309 


12 


26.3 


11 


7 


AO Men 


12 


38.5 


13 


9 


BD +2° 1729 


100 


9.8 


8 


4 


LQ Hya 


50 


18.3 


22 


8 


TWA 6A 


10 


77.0 


52 


25 


BD 1° 2447 


100 


7.2 


6 


3 


TWA 5B 


10 


50.0 


10 


7 


TWA 8A 


10 


21.0 


7 


3 


TWA 8B 


10 


21.0 


7 


2 


TWA 9A 


10 


64.0 


38 


19 


TWA 9B 


10 


64.0 


22 


13 


SAO 252852 


100 


16.4 


19 


10 


V343 Nor 


12 


39.8 


24 


12 


V2306 Oph 


100 


4.4 


4 


1 


HD 155555 AB 


12 


31.4 


22 


11 


HD 155555 C 


12 


31.4 


13 


8 


CD -64° 1208 


12 


29.2 


15 


9 


PZ Tel 


12 


49.6 


26 


15 


GJ 799 A 


12 


10.2 


4 


2 


GJ 799 B 


12 


10.2 


6 


2 


GJ 803 (Au Mic) 
BD -17° 6128 


12 


9.9 


4 


2 


12 


47.7 


25 


14 


GJ 813 


200 


13.6 


8 


4 


GJ 890 


100 


21.8 


26 


14 


HD 221503 


> 100-200 


13.9 


17 


9 



Table 5 
Statistical Results^) 



Cases AU Cases AU 



5Mj 


IOMj 


10 cases 


4 < r < 10 


9 cases 


1 < 


r 


< 5 


8 cases 


10 < r < 20 


8 cases 


5 < 


r 


< 10 


8 cases 


20 < r < 40 


9 cases 


10 < 


r 


< 20 


2 cases 


40 < r < 65 


2 cases 


20 < 


r 


< 36 



^ In the first and third columns the number of cases 
in which a 5Mj and a 10 Mj cxo-planct would have 
been detected in our survey. 
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